s.  ,  ,v  ■  '  . 


AD  A139057 


NUSC  Technical  Report  6883 
11  January  1984 


New  Formulas  That  Extend  Norton’s  Farfield 
Elementary  Dipole  Equations  to  the 
Quasi-Nearfield  Range 


Peter  R.  Bannister 

Submarine  Electromagnetic  Systems  Department 


4 

Naval  Underwater  Systems  Center  DTI  C 

Newport,  Rhode  Island  /New  London,  Connecticut  Ae’  Ills 

^  MAR  1  3  1384  S 

E 


§4  03  16  066 


Approved  for  pufcNc  release;  cUstributkm  cnNmHed. 


Preface 


This  report  was  prepared  under  N'JSC  Project  No.  A59007,  “ELF 
Propagation  RDT&E”  (U),  Principal  Investigator,  P.  R.  Bannister  (Code  3411), 
Navy  Program  Element  No.  1 1401 N  and  Project  No.  X0792-SB,  Naval  Electronic 
Systems  Command  Communications  Systems  Project  Office,  D.  Dyson  (Code  PME 
1 10),  Program  Manager  ELF  Communications,  Dr.  B.  Kruger  (Code  PME  1 10-XI). 

The  analysis  and  write  up  of  this  report  was  performed  while  the  author  was 
occupying  the  Research  Chair  in  Applied  Physics  at  the  Naval  Postgraduate  School, 
Monterey,  CA.  The  author  would  especially  like  to  thank  Professors  Otto  Heinz 
and  John  Dyer  and  Dean  Bill  Tolies  for  recommending  him  to  occupy  this  post  and 
NAVSEA  (Code  63R)  for  sponsoring  the  Chair. 

The  Technical  Reviewer  for  this  report  was  Anthony  Bruno. 


Reviewed  and  Approved:  1 1  January  1984 


(5 &  '0&eCdt_ 

D.  F.  Dence 

Head,  Submarine  Electromagnetic 
Systems  Department 


The  author  of  this  report  is  located  at  the 
New  Loudon  Laboratory,  Naval  Underwater  Systems  Center, 
New  London,  Connecticut  06320. 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTION'S 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMMK  2.  GOVT  ACCESSION  NO. 

TR  6883  AD*Ai3?  0*1  . 

1.  RECIPIENT* S  CATALOG  HUMSER 

4.  TITLE  umd  Smbtidti 

NEW  FORMULAS  THAT  EXTEND  NORTON'S  FARFIELD 

ELEMENTARY  DIPOLE  EQUATIONS  TO  THE 

QUASI -NEARFIELD  RANGE 

S.  TYPE  OF  REPORT  S  FERIOO  COVERED 

4.  PESPORMPW  ORG.  REPORT  NUMBER 

7.  AUTHOR** 

Peter  R.  Bannister 

L  CONTRACT  OR  ORANT  NUMBER*! 

1.  PERPOAMNG  ORGANIZATION  NAME  AMO  AOOAESS 

Naval  Underwater  Systems  Center 

New  London  Laboratory 

New  London,  Connecticut  06320 

10.  PROGRAM  CLEMENT.  PROJECT.  TASK 

AREA  SWORN  UNIT  NU  USERS 

A59007 

11.  CONTROLLING  OFFICE  NAME  AMO  A00RESS 

12.  REPORT  OATS 

11  January  1984 

13.  NUMREROPPAQES 

14.  MONITORING  AGENCY  NAME  S  A0ONCSS  lif  difftrmmt  frmmt  ComtrmUimf  Officmi 

IS.  SECURITY  CLASS.  M/iAm  rtpmrn 

UNCLASSIFIED 

ISA  OCCLASSMKATIONIOOWNORAMNQ 
SCHEDULE 

II.  OISTRItUTtON  STATEMENT  M/rt u  Htp*nt 


Approved  for  public  release;  distribution  unlimited. 


17.  OISTRISUTION  STATEMENT  tmf  thm  akmmct  tatarmd  im  Block  M.  if  diffmrmm  fnmm  Kapmrtt 


IS.  SUPPLEMENTARY  NOTES 


IE  KEY  WOROS  iCmamimmm  mm  rmm  udm  if  nmnmnry  rnrnd  idmmaify  by  kimck  mmmbmn 

Air-to-Air  Propagation  Horizontal  Magnetic  Dipole 

Electromagnetic  Fields  Vertical  Electric  Dipole 

Horizontal  Electric  Dipole  Vertical  Magnetic  Dipole 


21  AtSTftACT  tCvntimma  m  mum  if  mrwiiry  md  identify  by  btmek  mtmkmt 

New  formulas  for  the  electric  and  magnetic  fields  produced  by  the  four 
elementary  dipole  antennas  have  been  derived  for  the  air-to-air,  subsurface- 
to  air,  air-to-subsurface,  and  surface-to-surface  propagation  cases.  These 
formulas  are  of  rather  simple  form  and  reduce  to  previously  derived  results 
when  either  (1)  the  measurement  distance  is  much  less  than  a  free-space 
wavelength,  (2)  the  Soramerfeld  numerical  distance  is  small,  or  (3)  the 
measurement  distance  is  much  greater  than  a  free-space  wavelength.  Thej rjb&s-i. 
are  valid  at  any  frequency  and  at  any  range  beyond  a  certain  minimum  distance  — _ 


00,2*1473 


.20.  (Cont'd) 


f°2|the,flat>earth  Case'  1116  main  restrictions  on  these  formulas  are  (1) 
ln  I  i  (2)  the  measurement  distance  is  >.10  skin  depths  from  the  source, 
and  (3)  the  measurement  distance  is  >5  times  the  depth  of  burial  of  the  trans¬ 
mitting  or  receiving  point  sources. 

In  terms  of  computer  time,  these  new  formulas  can  be  evaluated  in  fractions 
of  a  minute  compared  with  hours  for  the  complete  numerical  evaluation  of  the 
exact  Sommerfeld  integrals. 

These  formulas  are  intended  to  supplement  the  author's  recently  derived 
subsurface- to- subsurface,  subsurface-to-surface,  surface-to-subsurface,  and 
surface-to-surface  propagation  formulas. 


TR  6883 


TABLE  OF  CONTENTS 


Page 

LISP  OF  TABLES . ii 

GLOSSARY  OF  SYMBOLS . iii 

INTRODUCTION  .  1 

AIR-TO-AIR  PROPAGATION  DERIVATION  PROCEDURE  .  2 

SUBSURFACE-TO-AIR  PROPAGATION  .  9 

AIR-TO-SU BSURFACE  PROPAGATION . 10 

SURFACE-TO-SURFACE  PROPAGATION  .  10 

DISCUSSION . 10 

CONCLUSIONS . 12 

REFERENCES . 43 


TR  6883 


LIST  OF  TABLES 


Table  Page 

1  Electric-Field  Air-tj-Air  Propagation 

Formulas  (jn2|  >_  10,  R.  >_  106)  . 15 

2  Magnetic-Field  Air-to-Air  Propagation 

Formulas  (|n2|  ^10,  Rj  ^106)  . 17 

« 

3  Electric-Field  Air-to-Air  Propagation  Formulas  for 

the  Farfield  Case  (|n2|  >_  10,  5y0ri  !  >>  1) . .  .  21 

4  Magnetic-Field  Air-to-Air  Propagation  Formulas  for 

the  Farfield  Case  (|n2|  10,  jy0R|  >>  1) . 23 

5  Electric-Field  Subsurface-to-Air  Propagation 

Formulas  (|n2|  ^10,  R  ^  106,  R  >_  5h) . . . 25 

6  Magnetic-Field  Subsurface-to-Air  Propagation 

Formulas  (|n2|  ^10,  R  106,  R  >_  5h) . 27 

7  Subsurface-to-Air  Propagation  Formulas 

When  p  »  z  (|n2|  ^10,  p  >_  i06,  p  5h) . 29 

3.  Subsurface-to-Air  Propagation  Formulas 

When  z  »  p  (jn2|  10,  z  106,  z  >_  5h) . 31 

9  Electric-Field  Air-to-Subsurface  Propagation 

Formulas  (|n2|  >^10,  D  ^  106,  D  5z) . 33 

.10  Magnetic-Field  Air-to-Subsurface  Propagation 

Formulas  (|n2|  ^  10,  D  106,  D  >_  5z) . 35 

11  Air-to-Subsurface  Propagation  Formulas 

When  p  »  h  (|n2|  >_  10,  p  >_  106,  p  ■>  5z) . 37 

12  Air-to-Subsurface  Propagation  Formulas 

When  h  »  p  (|n2|  >10,  h  >_  106,  h  >_  5z) . 39 


13  Surface-to-Surface  Propagation  Formulas 
for  p  >_  105  (|n2|  ^10,  z  =  h  =  0+) 


41 


D 

En 


“♦ 

Ez 

F 

h 


BED 

HMD 


H„ 


H, 


<fr 


H, 


I 

j0aP) 


m 

n 

n. 


n. 


P 

P 

R 


S£<WW5C»ee^^ 


TR  6883 


GLOSSARY  OF  SYMBOLS 


sin 
sin 


V¥(w)  /1  +  r«*\  f1  -  r..\ 

&  V,  ■  ■•■,  -I— H— h 

sin2  ♦1  -  42F(w)  /i  *  r„\  A  -  r„\ 

'An*,  +  4,  =  \~T~7!in  ♦*  '  \-T-jsiFM  *  si"  *i  -  V 


(p2  +  h2)1/2  (meters) 

Horizontal  electric-field  component  in  the  p  direction  (volts/meter) 
Horizontal  electric-field  component  in  the  <j>  direction  (volts/meter) 
Vertical  electric-field  component  (volts/meter) 

F(w)  or  F(w0),  Sommerfeld  surface-wave  attenuation  factors 


Height  (h  0)  of  transmitting  antenna  with  respect  to  earth's  sur¬ 
face  (meters) 


Horizontal  electric  dipole 
Horizontal  magnetic  dipole 

Horizontal  magnetic- field  component  in  the  p  direction  (amperes/meter) 
Horizontal  magnetic- field  component  in  the  <j>  direction  (amperes/meter) 
Vertical  magnetic-field  component  (amperes/meter) 

Current  (amperes) 

Bessel  function  of  the  first  kind,  order  zero,  with  argument  Ap 
Magnetic-dipole  moment  (ampere-meters2) 

V 

Yi/Y0,  index  of  refraction 
120ir,  free  space  impedance  (ohms) 
n(l  -  A2  cos2  lpj)1/2 
Electric-current  moment  (ampere-meters) 

Wait  integral 

(p2  ♦  z2)1/2  (meters) 
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u„ 


u 


1 

VKD 


VMD 

to 

X 


Yo 


A 

A, 


X 

Ar 


XV 


[p2  +  (z  -  h)2]1/2  (meters) 

LP2  +  (z  «■  h)2]1/2  (meters) 

Expf-YgR^/Rj ,  Sommerfeld  integral 
Time  (seconds) 

(A2  +  Yq)1/2  (meters-1)  (air) 

(X2  +  y^)1/2  (meters-1)  (earth) 

Vertical  electric  dipole 
Vertical  magnetic  dipole 
w0  or  w',  Sommerfeld  numerical  distances 
1  -  A2  cos2 


Height  (z  j>  0)  of  receiving  antenna  with  respect  to  earth’s  sur¬ 
face  (meters) 


sin  if>x  -  Aj 


,  Fresnel  reflection  coefficient  for  vertical 
sin  4 !»-  +  Ax  polarization 


sin  ipj  -  n2 

— -  ,  Fresnel  reflection  coefficient  for  horizontal 

sin  +  nj  polarization 


(-^PqEq)1/2  =  i2ir/X0,  upper  half-space  (air)  propagation 
constant  (meters-1) 


(itotijOj  -  a)2y j) 1/2,  lower  half-space  (earth)  propagation 
constant  (meters-1) 


Yo/Yj  s  1/n 

A(1  -  A2  cos2  ij/j) x/2 


/  2  X1/'2  /to2e2  X1/2  oje1 

\“»*o aJ  L\  °i  /  °i. 


-1/2 


,  skin  depth  in  the  water  or 
earth  (meters) 


=  10-9/36tt  farads/meter,  permittivity  of  free  space 
Permittivity  of  lower  half-space  (earth)  (farads/meter) 

Rummy  integration  variable  in  the  basic  Sommerfeld  integrals  (meters-1) 
Wavelength  in  free  space  (meters) 
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p  ( x 2  +  y2)1/2  radial  distance  in  a  cylindrical  coordinate 

system  (meters) 

Oj  Conductivity  of  the  lower  half-space  (earth)  (Siemens/meter) 

<(>  tan-1  (y/x),  azimuth  angle  in  a  cylindrical  coordinate  system 

y  -  nQ  =  4ir  10~7  henries/meter,  permeability  of  free  space 
Lan-1  (z/p)  or  tan-1  (h/p),  elevation  angle 
i|»0  tan-1  [(z  -  h)/p],  elevation  angle 

ifij  tan-1  [(z  +  h)/p],  elevation  angle 

to  2irf  radians/ second,  angular  frequency 
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3  M 


'S  % 

2  1 


v/vi 

Reverse  Blank 


-  av^r ,  &&  r  . 


TO  6883 


NEW  FORMULAS  THAT  EXTEND  NORTON'S  FARFIELD  ELEMENTARY 
DIPOLE  EQUATIONS  TO  THE  QUASI -NEARFIELD  RANGE 

INTRODUCTION 


It  is  the  purpose  of  this  report  to  present  new  formulas  for  horizontal 
electric  dipole  (HED) ,  horizontal  magnetic  dipole  (HMD) ,  vertical  electric 
dipole  (VED) ,  and  vertical  magnetic  dipole  (VMD)  air-to-air,  subsurface-to- 
air,  air-to-subsurface,  and  surface-to-surface  propagation.  The  new  air-to- 
air  propagation  formulas  extend  Norton's  results,! >2,3  which  are  stated  to  be 
valid  for  measurement  distances  greater  than  a  free-space  wavelength  (AQ), 
down  to  the  quasi-nearfield  range,  which  is  defined  as  the  range  where  the 
measurement  distance  is  much  less  than  a  free-space  wavelength  but  much  greater 
than  an  earth-skin  depth  (6) .  The  new  subsurface-to-air  and  air-to-subsurface 
propagation  formulas  reduce  to  previously  derived  results^-lO  when  the  meas¬ 
urement  distance  is  much  less  than,  or  comparable  to,  XQ.  Norton's1 >2,3  (cor¬ 
rected)  air-to-air  propagation  formulas  are  summarized  by  King,11  while 
Kraichman12  has  tabulated  the  subsurface-to-air  and  air-to-subsurface  propa¬ 
gation  equations  derived  by  various  authors. 4-10 

In  the  past,  many  investigators  erroneously  have  believed  that  the  field- 
strength  equations  tabulated  in  Chapter  3  of  Kraichman12  are  only  valid  when 
the  conduction  currents  in  the  water  or  earth  are  much  greater  than  the  dis¬ 
placement  currents  (i.e.,  »  wej).  Indeed,  as  long  as  |n2|  =  | Yf /y© i  >:>  !> 

the  displacement  currents  can  be  included  simply  by  replacing  exj  by  o1  +  iuxej 
in  the  field-strength  equations.  Thus ,  Kraichman' s  tabulated  results  are  con¬ 
siderably  more  general  than  they  are  stated  to  be. 

The  formulas  presented  in  this  report  are  intended  to  supplement  the 
author's  recently  derived  subsurf ace-to-subsurf ace,  subsurface-to-surface, 
surface-to-subsurface,  and  surface-to-surface  propagation  formulas. !3  They 
are  valid  at  any  frequency  and  any  range  beyond  a  certain  minimum  distance  for 
the  flat-earth  case.  The  main  restrictions  on  these  formulas  f»'.e  (1)  |n2|  >^ 

10,  (2)  the  measurement  distance  is  >106  from  the  source,  and  (3)  the  measure¬ 
ment  distance  is  >5  times  the  depth  of  burial  of  the  transmitting  or  receiving 
point  sources. 

For  the  air-to-air  propagation  case,  the  four  dipole  antennas  (VED,  VMD, 
HED,  and  HMD)  are  situated  at  height  h  (h  >_  0)  with  respect  to  a  cylindrical 
coordinate  system  (p ,4>, z)  and  are  assumed  to  carry  a  constant  current,  I.  The 
axes  of  the  VED  and  HED  (of  dipole  moment  p)  are  oriented  in  the  z  and  x 
directions,  respectively,  while  the  axes  of  the  VMD  and  HMD  (of  dipole  moment 
m)  are  oriented  in  the  z  and  y  directions,  respectively.  The  earth,  which  is 
assumed  to  be  a  homogeneous  medium  with  conductivity  o1  and  dielectric  con¬ 
stant  £j  (=  er~  ) ,  occupies  the  lower  half-space  (z  <  0)  and  the  air  occupies 
the  upper  half-space  (z  >  0) .  The  magnetic  permeability  of  the  earth  is 
assumed  to  equal  yQ,  the  permeability  of  free  space.  Meter-kilogram-second 
(MKS)  units  are  employed  and  a  suppressed  time  factor  of  exp(iwt)  is  assumed. 
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AIR-TO-AIR  PROPAGATION  DERIVATION  PROCEDURE 


As  an  example  of  our  derivation  procedure,  consider  an  HED  source.  When 
h  and  z  are  >0,  the  Sommerfeld  integral  expressions  for  the  HED  Hertz  vector 
are4»5,14 


nx  = 


p . 

4iriwe0 

Ro 

R1 

Ri  Jn  \ui  +  uo/ 


0(z+h)  i  1 

J0(Xp)-<iX 


(1) 


and 


p  cos  *  3  f  2(ul  05  “u0(z+h)T  „ 

n  =  - SL  X  —  I  -5 - 5 —  e  u  JQ(Xp)XdX  , 

4irxa)e0  3p J  YfuQ  +  Y^i  0 


(2) 


where 


Rg  =  p2  +  (z  -  h)2. 


R2  =  p2  +  (z  +  h)2, 
u0  =  xZ  +  Y0’ 


u2  =  X2  4-  y2, 

Yo  =  '“2P0e0«  311(1 

Y2  *  +  iuEj). 

From  equations  (1)  and  (2) ,  and  utilizing  the  identity  (Uj  -  uQ) (uJ  +  uQ) 


p  cos  4> 

3 

e'raRo 

e"YoRi 

4iriu)e0 

‘  8P 

.  Ro 

R1 

.»2y2e-VZ+h) 


°jl +r_!£__ 
1  Jn  ylU0  + 


Y0U1 


Jn(Xp)XdX 


(3) 


When  |n2|  »  1  and  ReCYjRj)  »  1,  the  function  u:  in  the  exact  integral 
expressions  can  be  replaced  by  Yj*  the  propagation  constant  in  the  earth. 4 
Therefore, 


"*  - 


P 

VY»R» 

4iriue0 

.  R° 

‘Y0Rl  2  T  -u.(z+h) 

-r-  *  TTTZkJ  (yr-  0  JoCxo)“« 

Ri  tYi  r0JJn 


(4) 
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Since  Sommerf eld's  integral,  Sj,  is  equal  to 


o 

A 


-un(z+h)  x 

e  0  J0(Xp^X 


(T1 


(5) 


then, 


A- 


n (z+h)  3Sl 

0  Jn(Ap)AdA  =  - 


S(z  +  h) 


si"  -y.r, 

—j-tl  *  Y0R,)e  «  > 


(6) 


and 


oo 

/" 


-un(z+h)  a2si 

,e  0  J0(Xp)XdX  = 


3(z  +  h)2 


(7) 


r3i 


[(1  +  YqRj)  -  Sin2  ^(3  +  SYqRj  +  YjRf)]  > 


where  sin  <fix  =  (z  +  h)/Rj.  Therefore,  from  equations  (4),  (6),  and  (7), 

-YnR, 


"x  ' 


<!  =  (Z 

+  h)/R 

P 

e_Y°R° 

4iriwe0 

Ro 

2e 


'0”1 


R1  (T?  -  Y§)R? 

[(1  +  y  1^1  sin  <^3(1  +  y0Rx)  -  sin2  <1^(3  +  3y0Rx  +  Y§Rf)]  J 


(8a) 


Alternatively,  we  could  have  replaced  the  quantity  Uj  in  the  exact  inte¬ 
gral  expressions  by  yx[1  -  (Yq/yx)cos2  ij^]1/2  instead  of  y^3’*1’*5  This 
change  is  equivalent  to  setting  the  normalized  surface  impedance,  Aj,  equal  to 

Ax  =  A(1  -  A2  cos2  ijij)1/2  ,  (8) 

where  A  =  Yq/Yj  =  1/n  and  cos  ^  =  p/Rx. 

Since  (y2|  »  |yq|»  the  factor  in  parentheses  in  equation  (9)  is  very 
near  unity  (i.e.,  A}  ~  A) .  Thus,  the  resultant  modification  is  usually  small. 
The  retention  of  the  parenthesized  term  is  justified  by  the  fact  that  the 
exact  form  of  the  Fresnel  reflection  coefficient  is  recovered  from  the  asymp¬ 
totic  solution, For  example,  if  we  utilize  the  identity  (ux  -  u0)(Uj  +  u0) 
=  Y2  “  Yq>  substituting  ux  =  yx(1  -  A2  cos2  i^)1/2  into  equation  (1)  and  eval¬ 
uating  the  Sommerfeld  integrals  results  in 
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Y 

w  =  -  -~{sin  1|<1  +  A)2  . 


If  we  replace  A  by  (equation  (9)), 


/kl\1/2  f1  '  r,l\  =  Ai 

\w/  \  2  /  sin  +  Aj  * 


where 


sin  *  -  A 

r  =  - - i - i-  (2£ 

sin  ij»1  + 

is  the  Fresnel  reflection  coefficient  for  vertical  polarization.  Therefore, 
equation  (21)  can  be  expressed  as 


where 


[1  -  F(w )]'- 


F(w)  -  1  -  i(irw)1/2e“werfc(iw1/2) 


is  the  Sommerfeld  surface-wave  attenuation  function  and 

4Wt 

w  =  _  JL  (sin  ^  ♦  Ax)2  =  —  —^2 

is  the  Sommerfeld  numerical  distance.  For  small  numerical  distances  F(w)  -  1 
and  for  large  numerical  distances  and  negative  arguments  F(w)  ~  -l/(2w). 


Since 


sin  ^  +  AjF(w) 
sin  ljij  +  L1 


=  A  , 


equation  (20)  reduces  to 
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Another  factor  that  we  will  encounter  in  the  derivation  of  the  field- 
strength  components  is  the  factor  B,  which  is  equal  to 


B  =  sin  ^  -  AjA 


(i  *  r„\  /i  -  r„\ 

=  — ^sin  -  ^-y-jAiFCw) 


sin2  iJ/j  -  A2F(w) 
sin  ^  *  &1 

For  small  numerical  distances  (i.e.,  F(w)  -  1),  A  ~  1,  and  B  -  sin  ^ 

-  Aj.  Furthermore,  for  sin  »  JaJ,  A  ~  1  and  B  ~  sin  •  For  sin 
comparable  to  or  less  than  Aj,  the  horizontal  distance  p  will  be  much  greater 
than  the  sum  of  the  transmitting  and  receiving  heights  (z  +  h) .  In  the  limit 
as  ^  approaches  0,  A  ~  F(w0)  and  B  -  -AjFCWq),  where 

F(wfl)  -  1  -  i(irw0)1/2e  W°erfc(iw01/2)  (32) 


V  a2 

w0  “  ‘  2  1  ' 

When  p2  »  (z  +  h)2,  Wait4flS  has  shown  that 

F(w)  ~  [1  +  y0Mz  +  h)]F(w0)  =  [1  +  YoR^!  sin  4»1]FCw0) 


3F(w) 


ToaifCw0) 


Extensive  numerical  results  for  the  function  F(wQ)  have  been  provided  by 
Wait^S  and  King  and  Schlak.16 

Since  the  factor  A  (equation  (29))  is  different  from  unity  only  when  (1) 
the  angle  $  is  very  small  and  (2)  the  Sommerfeld  attenuation  function  F(w)  is 
different  from  unity,  A  is  only  a  farfield  surface-wave  term.  Therefore,  we 
can  discard  all  derivatives  of  A  that  are  not  farfield  terms.  For  example, 
when  |yf |  »  1  y§  I  311(1  |aJ  sin  ^  «  1, 


e“Y0Rl 

=  -  A(1 

+  YqR^COS 

♦l 

Ri 

e'YoRl 

-  -  A(1 

+  YqRjDCOS 

R? 

e'Y0Ri 

“  -  (H 

■  YnR,A)cos 

♦i 

o2 
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Therefore,  from  equations  (30)  and  (36), 
n  cos  A  r*r><;  tli  p>  -Yori  r 


p  cos  <t>  cos  i^e  0  1  sin  *», 

®z  - - T~. - Tt -  (1  ♦  YnR.A) - (3  +  3yAR, 

2niu»e0Y1Ri  0  1  YjRj  0  1 

When  |n2|  »  1  and  Re(Y1R1)  »  1,  equation  (3)  reduces  to 


p  cos  4>  cos  ^e  u  1  sin  ^ 

- - (1  +  Y0R1A)  -  7S~(3  +  3Y°R1  *  Y0R1)  •  (3?) 


.  1  ~ 


3 

e~Y°Ro 

3p 

Ro 

2  re-uo(z+h) 

4  /  5 - -  J0(Xp)XdX 

n2J0  uo  +  V  0 


Since 


l 

I 

rH  j 

1  \ 

u0  +  Y0A 

1 

C  1 

°  1 
1 

\u0 

u0  +  To A/ 

then,  following  the  same  procedure  as  in  the  derivation  of  the  HED  nz  compo¬ 
nent  results  in 


4iriu)ef 


_3_[e2 

8p[  R 


e'Y°R0 

e‘YoRl 

+  2  4 

.  Ro 

R1 

e-T°R0 

e_YoRi 

2Ae” 

-YnRn  "YnRi 

wf  C1  *  Y0R0icos  *0  ~r2  -  U  ♦  Vl>“S  *1  ~"r2 


+  —a  +  Y0RiA)cos  V 
n* 


e'T°R‘1 

1  «!  r 


Since  we  have  now  derived  expressions  for  the  HED  Hertz  vector  (equations 
(8),  (37),  and  (39)),  the  fields  in  air  can  be  obtained  from 

.  2  =  -Y$  +  1$  X  it) 

.*  *  *  C40) 

H  =  iue0(v  x  n)  . 

By  following  the  same  procedure  as  outlined  above,  we  can  also  obtain 
suitable  expressions  for  the  HMD,  VED,  and  VMD  Hertz  vectors.  The  resulting 
HED,  HMD,  VED,  and  VMD  field-component  expressions  for  the  air-to-air 
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propagation  case  are  presented  in  tables  1  and  2.*  They  are  strictly  valid 
for  jn2 |  >>  1  and  ReCyjR^  »  1.  However,  for  most  cases,  the  requirement  that 
|n2|  10  and  Rj  106  is  sufficient.  When  IyqR^ |  <<  1»  they  reduce  to  Bannis¬ 

ter’s  quasi-near  range  results.?*10 

When  JyqRJ  »  1,  they  reduce  to  Norton’s  farfield  range  results1*2*3*11 
(for  jn2 |  »  1).  For  convenience,  the  farfield  expressions  ( I Y0Ri |  >>  !)  are 
listed  in  tables  3  and  4.  In  these  tables,  sin  \pQ  =  (z  -  1i)/Rq,  cos  t|)0  =  p/Ro» 
sin  ^  =  (z  +  h)/Rj,  and  cos  4* j  =  p/Rj. 

The  air-to-air  propagation  results  presented  in  this  report  can  be  extended 
to  a  multilayered  earth  simply  be  substituting  y^Q  for  Yj*  where  Q  is  the 
familiar  plane-wave  correction  factor  employed  to  account  for  the  presence  of 
stratification  in  the  earth.15  For  a  homogeneous  ground,  the  argument  of  the 
Sommerfeld  numerical  distance  wQ  (equation  (33))  is  always  between  0  and  -90 
deg,  resulting  in  the  transverse  magnetic  (TM)  surface-wave  fields  varying  as 
1/p2  as  p  ■+<*>.  For  a  stratified  ground,  the  argument  of  wQ  can  be  positive, 
resulting  in  the  TM  surface-wave  fields  varying  as  1/Vp*  This  fact  as  dis¬ 
cussed  in  further  detail  by  Wait15  and  King  and  Schlak.10 


SUBSURF ACE-TO-AIR  PROPAGATION 


The  HED,  HMD,  VED,  and  VMD  field-component  expressions  for  the  subsurf ace- 
to-air  propagation  case  (h  £  0,  z  >_  0)  can  be  obtained  from  the  air-to-air 
propagation  equations  (tables  1  and  2)  simply  by  setting  h  =  0  and  multiplying 
each  resulting  expression  by  exp(Yjh).  (All  VED  components  must  also  be  multi¬ 
plied  by  1/n2  to  satisfy  the  boundary  conditions.)  The  resulting  equations  are 
presented  in  tables  5  and  6  for  the  general  case,  in  table  7  for  p  »  z,  and  in 
table  8  for  z  »  p.  In  these  tables,  R  =  p2  +  z2,  sin  i|>  =  z/R,  and  cos  <Ji  = 
p/R. 

The  expressions  presented  in  tables  5  through  8  are  strictly  valid  for 
|n2!  »  1,  Re(YjR)  »  1,  and  R  »  |h|.  However,  for  most  cases,  the  require¬ 
ment  that  | n2 1  10,  R  106,  and  R  >_  5|h|  is  sufficient. 

When  |yqR|  <<:  !>  they  reduce  to  Bannister's  quasi-near  range  formulas,?*10 
which  are  tabulated  in  Kraichman1^  (with  a1  replaced  by  Oj  +  icoe^.  When  F(w) 

=  1  (i.e.,  small  numerical  distances),  they  reduce  to  Bannister's  nearfield 
range  formulas,9*10  which  are  also  tabulated  in  Kraichman12  (with  Oj  replaced 
by  Oj  +  iwej,).  Furthermore,  when  h  =  0  and  |yoR|  >>  1»  they  are  consistent 
with  Norton’s  farfield  surface-to-air  results. i»2, 3, 11 


*A11  tables  have  been  placed  together  at  the  end  of  this  report. 
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AIR-TO-SUBSURFACE  PROPAGATION 


Vr. 


h 


The  HED,  1.MD,  VED,  and  VMD  field-component  expressions  for  the  air-to- 
subsurface  propagation  case  (h  ^  0,  z  <_  0)  can  be  obtained  from  the  air-to-air 
propagation  equations  (tables  1  and  2)  simply  by  setting  z  =  0  and  multiplying 
each  resulting  expression  by  expfyjZ) .  (All  Ez  components  must  also  be  multi¬ 
plied  by  1/n2  to  satisfy  the  boundary  conditions.)  The  resulting  equations 
are  presented  in  tables  9  and  10  for  the  general  case,  in  table  11  for  p  »  h, 
and  in  table  12  for  h  »  p.  In  these  tables,  D  =  p2  +  h2,  sin  \ (i  =  h/D,  and 
cos  ip  =  p/D. 

The  expressions  presented  in  tables  9  through  12  are  strictly  valid  for 
|n2 1  »  1,  Re(y1D)  »  1,  and  D  »  |z|.  However,  for  most  cases,  the  require¬ 
ment  that  |n2j  >_  10,  D  106,  and  D  5|z|  is  sufficient. 

7  10 

When  |y0D|  «  1,  they  reduce  to  Bannister's  quasi-near  range  formulas,  * 
which  are  tabulated  in  Kraichmanl2  (with  al  replaced  by  Oj  +  itocj).  When  the 
numerical  distance  is  small  (i.e.,  F(w)  ~  1),  they  reduce  to  Bannister's  near¬ 
field  range  formulas, 8  which  are  also  tabulated  in  Kraichmanl2  (with  Oj 
replaced  by  Oj  +  iajex) .  Furthermore,  when  z  =  0  and  |y0D|  »  1,  they  are  con¬ 
sistent  with  Norton's  farfield  air-to-surface  results. 1,2, 3, 11 


SURFACE-TO-SURFACE  PROPAGATION 


The  simple-form  HED,  HMD,  VED,  and  VMD  field-component  expressions  for 
the  surface-to-surface  propagation  case  can  be  obtained  from  the  air-to-air 
propagation  equations  (tables  1  and  2)  simply  by  setting  both  z  and  h  equal  to 
zero.  The  resulting  equations  are  listed  in  table  13.  They  are  strictly  valid 
for  jn2|  »  1  and  Re(y1p)  »  1.  However,  the  requirement  that  |n2|  >.  10  and 
p  ^106  is  sufficient.  When  either  (1)  Jy0p  |  «  1,  (2)  F(w0)  ~  1,  or  (3) 

|y0p|  »  l,  they  reduce  to  previously  derived  results. 4-12  Note  that  the 
function  F  in  this  table  is  equal  to  F(wQ). 


DISCUSSION 


Since  |n2|  »  1  (i.e.,  | A2 |  «  1),  the  factor  nx  -  n  and  Ax  ~  A.  Fur¬ 
thermore  , 

sin  -  u 

r„ - - -  (4D 

sin  +  A 


and 

sin  ij),  -  n  2  sin  ^ 

r  .  - - 1  +  - 

sm  ^  +  n  n 

10 


(42) 


CffUfri  '•X'U*  'ttr+'Sl&m  *&$&&V2S%ZgN&rtf*p  *-v-  +*&■* 
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For  small  numerical  distances  (i.e.,  F(w)  -  1) ,  A  ~  1  and  B  ~  sin  \p^  -  A. 
Furthermore,  for  sin  tp^  »  |a|,  A  ~  1  and  B  -  sin  <pj.  When  sin  is  compa¬ 
rable  to  or  less  than  A,  the  horizontal  distance  p  will  be  much  greater  than 
the  sum  of  the  transmitting  and  receiving  antenna  heights  (z  +  h) .  In  the 
limit  as  ipj  approaches  zero,  A  -  F(wfl)  and  B  -  -AF (wQ) ,  where 


F(wn)  -  1  -  i(Trw0)1/2e"W°erfc(iwJ/2) 


V  a2 

w0  =  "  —  A2  * 


For  this  case  (i.e.,  p2  »  (z  +  h)2).  Wait4*15  has  shown  that  F(w)  can  be 
replaced  by 


F(w)  -  [1  +  y0A(z  +  h)jF(wQ)  , 


and  we  can  make  use  of  his  tabulated  resultslS  of  the  function  F(wQ) 


When  we  were  confirming  the  validity  for  the  subsurface-to-air  and  air- 
to-subsurface  propagation  equations  derived  in  this  report  (which  were  derived 
from  the  air-to-air  propagation  equations) ,  we  discovered  that  they  could  also 
have  been  obtained  by  combining  previously  derived  results.  The  derivation 
procedure  can  best  be  shown  by  example. 


From  table  3.3  of  Kraichman,12  the  nearfield  range  subsurface-to-air  HED 
component  is 


Yih  -Y0R 

Hf  -  -  ■»  *  ^  *  4»2>  • 


When  |y0R|  »  1>  equation  (46)  reduces  to 


Y,h  -yqR 

Y nP  cos  <{>e  A  e 

hHE  .  121 _ _ 

H<P  2ttR 


(A)  . 


If  we  take  Norton's  farfield  equation  (table  4),  set  h  =  0,  and  multi¬ 
ply  by  exp(Yxh),  we  obtain  (for  |y0R|  >>  1) 


Hf  -- 


Y0P  cos  <pe 


Tih  -YftR 


.n  -YnK  /  \ 


[sin  \p  +  AF(w)] 


Since 


l1  “  F"\  A 

\  2  /  sin  ip  +  A 
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then. 


(— j— )[sin  *  ♦  AF(wJJ  .  JSa.t  t  *■*»>.] 
\  2  /  l  sin  ip  *  A  J 


Therefore,  equation  (48)  reduces  to 


Y„h  -y„R 
1 1  o  0 


Ynp  cos  <f>e  1  e 


— (AA)  , 


which  is  equivalent  to  equation  (47)  except  for  the  factor  A.  When  either 
sin  tj/j  »  J  A  j  or  the  Sommerfeld  numerical  distance  (w)  is  small,  A  -  1.  For 
|n2|  »  1,  the  range  of  validity  of  equations  (46)  and  (51)  overlap  when  |w[ 

«  1  and  |y0R)  »  1  simultaneously.  Therefore,  we  can  simply  combine  equations 
(46)  and  (51)  to  obtain  an  expression  for  the  HED  component  valid  from  the 
quasi-nearfield  to  the  farfield  ranges.  Therefore,  for  |n2|  »  1,  Re(yiR)  » 

1,  and  R  »  |h } , 


jiK'roR 


2TTYJR3 


(1  +  YftR  +  Y^R2A)  , 


which  is  the  result  listed  in  table  6. 


In  following  this  procedure,  we  discovered  a  small  error  in  the  HED  and 
HMD  E^  and  Hp  farfield  subsurface-to-air  and  subsurface-to-subsurface  propa¬ 
gation  equations  tabulated  in  tables  3.1  and  3.2  of  Kraichmanl2  (and  also  in 
equation  (14)  of  Wait4) .  These  equations  should  be  multiplied  by  the  factor 
[1  +  F(wq)]/2.  They  will  then  agree  with  Norton's  farfield  surface-to-surface 
propagation  results.  This  correction  is  unimportant  for  small  numerical  dis¬ 
tances,  since,  for  this  case,  F(w3)  -  1  and  [1  +  F(wQ)]/2  =  1.  For  large 
numerical  distances,  the  difference  is  a  factor  of  1/2.  however,  for  this 
case,  Ep  »  E^  and  »  Hp. 


CONCLUSIONS 


New  formulas  that  extend  Norton's  farfield  elementary  dipole  equations  to 
the  quasi-nearfield  range  have  been  developed  for  the  air-to-air,  subsurface- 
to-air,  air-to-subsurface,  and  surface-to-surface  propagation  cases.  They  are 
valid  at  any  frequency  and  at  any  range  beyond  a  certain  minimum  distance  for 
the  flat-earth  case.  The  main  restrictions  on  these  formulas  are  (11  |n2|  >_ 
10,  (2)  the  measurement  distance  is  >10  skin  depths  from  the  source,  and  (3) 
the  measurement  distance  is  >5  times  the  depth  of  burial  of  the  transmitting 
or  receiving  point  sources. 

These  new  formulas  reduce  to  previously  derived  results  when  either  (1) 
the  measurement  distance  is  much  less  than  a  free-space  wavelength,  (2)  the 
Sommerfeld  numerical  distance  is  small,  or  (3)  the  measurement  distance  is 
much  greater  than  a  free-space  wavelength. 
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It  should  be  roted  that  the  two  media  can  be  inverted  and  the  air 
replaced  by  the  earth's  crust  (of  conductivity  o2  and  dielectric  constant  e2) . 
The  same  equations  (tables  1  through  13)  can  be  utilized  as  long  as  |n|j  = 
|y2/y| |  >.10,  Rj  >_  1062*  and  R  >_  5|h|  (or  D  >.  S|z|)  simply  by  replacing  iueQ 
by  o2  *  i(oe2« 

We  have  recently  employed  finitely  conducting  earth-image  theory  tech¬ 
niques  to  derive  field-component  expressions  valid  at  any  range  from  the  source 
for  the  air-to-air  propagation  case.  The  only  restriction  on  these  new  image- 
theory  formulas  is  that  VI  »  1.  When  the  measurement  distance  is  >10  skin 
depths  from  the  source,  they  reduce  to  the  air-to-air  equations  derived  in  this 
report.  These  new  image-theory  results  will  be  the  subject  of  a  future 

report. 17 
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Table  1.  Electric-Field  Air-to-Air  Propagation  Fd 


Dipole 

Type 


4*iue„ 


(3  ♦  3y0R0  +  Y§R§)sin  *Q  cos  *0  “p“ 

Ko 


♦  (3  ♦  3y0R1  ♦  r„YjR2)sin  +:  cos  *x  —p 

2y,  cos  *je~v°Rl  f  /l  -  r„\  11 

-  -  :  --n^  "I1  +\— )F(w)YoRiJJ 


4*iwer 


Vf"S~~*lt(3  cos2  *0  -  D(1  ♦  Y0R„)  -  Y§RJi  Sin2  *0J— 


0  O'  '(TO  r0-‘  d3 

0 

e”Y0Rl 

_  OnO  .9  «t© 


-  [(3  cos2  ^  -  1)(1  +  YqRj)  -  r„YgRf  Sin2  *11~3  — 
+  2"724~'[1  ‘  YiRi  sin  *i  +  YoRi  +  ("~7“)f(w)yori]} 


i»u0*  cos  4  j  -Y0R0 

_ i _ If*  x  n  j.  c 


iuy0“  [  e"Y°R° 


-  p  +  y0r0)cos  *o  —$r-  -  f1  i 

-y.r, 

2  cos  ♦.e  0  1 

+  (Y?  -  Yq)Ri  +  YlRl  *in  *1^3  ^ 

-  sin2  (15  ♦  ISyqR!  +  6YqR2  +  Y^Rf)l 


_  J(1  +  Y0R0)sin  +0  ^ 

e‘YoRi 

.♦  <1  4  VoR^sin  5_|— 


i«P0®  sin  *  f  "yoRo 

__!(l ,  Y/o)5in  ..y, 

♦  ^r^r-[2  *  W1  *  *)  -  sin*  *,(j 
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'vww3w»n* » ’j^Bsggjw 


Air  Propagation  Formulas  (|n2j  >_  10,  Rj  >_  106) 


t(l  '  3  Sin2  W  +  Y0R0)  ♦  Y2R2  COS2  *0]~3- 


+  [(1  -  3  sin2  Cl  +  YqR1)  +  r(tYoRl  cos2  *il“ 

~^0R1 

+  (1  -  r„)F(w)Y2R2  cos2  1|«1  — ^| — 


|  -YoR0  ”^ori 

Is  +o  -  u +  y0ri)cos  *i  ^-p— 

i;  Ko  Ri 


I*  Y^  sin  (3  ♦  SYjjRj  +  YqRj) 


||  +  6Y2Rf  +  YgRf)] 


W-p-  -  a  ♦  Y0Ri  ♦  Y§Rf)^-p 

*0  Kj 


|sin  ♦j)  +  y0Rj  (1  +  A  +  y^  sin  i^) 


Y§Rf)l 


-£■  “((3  *  **  *  T«fi")sin  *0  “s  *0  ^ 


(3  +  SyqRj  ♦  r„Y§Rf)sin  cos  ^  2_— 
2Tj  cos  *;=  T»Rlf  fl  -  r„^  ]| 

— Sf — [‘  *  (— >mv.]J 


iuu0m  cos  $ 


,  .  „  ,-V.  ."'oRi 

>sin  *0  -p~  *  C1  *  Y0Rl)Sin  *1  V" 
0  R1 

A)  -  sin2  ^(3  ♦  Sy^j  +  YqR|)3 


(1  ♦  Y0Ro)cos  *„ 


♦  (1  ♦  rMYoRi)cos 


,“yori 


e"YoRi 

♦  (1  -  r(l)F(w)Y0Rj  cos  *j  — ^5— 
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Table  2.  Magnetic-Field  Air-to-Air  Propagation 


Dipole 

Type 

«P  ‘  . 

«♦  ! 

VED 

0 

p  f  e”Y°R° 

iju  .  ,0R0)COS  ♦„  .(It 

e"YoRij 

4  (1  -  r„)F(w)Y0R1  cos  +i  ‘  Rz  1 

4 

f 

VMD 

B  e"Y°R° 

-|(3  ♦  3YoR0  ♦  Y2R2)sin  +o  cos  *o  R3 

“Yo^l 

-  (3  ♦  3y0R,  +  YoRf)s^n  *x  cos  ^  r3 

2  cos  f.e  Y°Rl 

K3  +  3YoRl  +  YoRi) 

-  sin2  ♦jCIS  ♦  15y0R1  +  6y2R2  +  YqRj)] 

% 

'i 

4 

i 

§ 

% 

1 

•% 

0  1 

1 

=?. 

% 

1 

i 

ft 

i 

c lMNM«4>i ; «a*M \+,  steft  lif  .tho*!**  1 


Dipole 

Type 


(Cont'd)  Magnetic-Field  Air-to-Air  Pro 


-  P  (1  ♦  Y0R0^sin  *0 


e"Y0R0  e'Y°Rl 

— -  (1  +  YoR^sin  i>l  — jjj— 


?*‘Y0R1 

y~r| — [2  +  Y°Rl(1  +  A)  ■  sin2  4(1(3  +  3y°Ri  *  YoRi}1 


y0R0)sm  *0 


(1  +  ^.YoR^sin  — -r 


!-^LA  [2d  ♦  Y0r0)  -  sin2  *0(3  ♦  3y0R0  ♦  Y^)]^ 

-yori 

♦  [2(1  ♦  YoRjA)  -  sin2  ^(3  ♦  SYoRj  +  Y^Rf)]— — 


-  1(1  ♦  Y0R0  ♦  Y^)5- 

-y0R1 

+  (1  +  YqRi  +  ri|YoRi)  p 
,  ,  e'YoRil 

+  (1  -  r„)F(w)Y2R2  -^3-j 
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Dipole 

Type 


Table  3.  Electric-Field  Air-to-Air  Propagation  Formulas  for  the-; 


w . 


,-y0ro 


— |sin  ipQ  cos  *0  — —  +  r„  sin  ^  cos  ^  — 


-  (1  -  r„)AF(w)cos  *  5_ 

R1 


Y0n0p  cos  *  [  e-Y0R0 

- n - sinZ  —7 — 


o  nr-  -  r« sin2  *i  V 
*0  K1 


-Y0R1 


-  (1  -  r„)A2F(w)- 


-vr 

R1  . 


iojp0Yom  cos  *  [  e"Y°R°  e' 


s m  *0  —5—  ♦  r„  sin  *t 


(i  -  r„)AF(w)- 


r*oRl 


x“>Vo 


— [COS  *0  —  -  cos 


2  cos 

Y0nZRi 


Y0n0p  sin  * 


2e‘Y°Rl  1 

*  *  *  *  TA  «  *l)j 


icoiioYo®  sin  ♦  e~Y°Ro  1 

- -  sin  ilin  - + 1 

Att  V°  R0  I 


-^°Ri r  ,  1 

1  1  y0R 

» 
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agation  Formulas  for  the  Farfield  Case  (|n2|  10,  |yqR1 j  »  1) 


P  sin  *Je'YoRo  e'Y0Rl 


L-Y0ri 


(1  +  A  +  YjRj  sin  +!) 


Y0nQp  cos  t  -Y0r0 

- : - sin  i>„  cos  *„  - 

4ir  Rg 

e"YoRl 

-  r„  sin  *j  cos  ^  — — 


+  (1  -  r„)AF(w)cos  1>1 


»"Y0R1 


|*0®  sin  ♦ [  e-YoRo  e“Y0Rl 


n; — -+in  *0  — ♦ 


i«u0Y0“  cos  +  e-Y0R0  -YqRj 

— * — rs  *°  ~ + r”  cos  * — 


boRlf.,  ! 

NrH  *>  ■  v 


Cl  +  A  -  3  sin2  *j)  II  +  (1  -  r,|)F(w)cos  4> 


e~Y°Ri 

1  Ri  . 
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Propagation  Foraulas  for  the  Farfield  Case  (|n2|  >_  10,  J yqR|  »  1) 
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Dipole 

Type 


II 

II 


Table  S.  Electric-Field  Subsurf ace-to-Air  Propagation  Forau 


-Y,h  -Y0R 
cos  *e  1  e  0 


2*(oj  +  iwe^R3 


[(3  +  3y0R)sin  * 


Ay„R  ♦  YqR2®] 


cos  *e 


-Y,h  -Y„R 


2*(o1  +  iuejJR3 
►  Y0r  *  Y^RznB) 


[1  -  YjR  sin  * 


-y  h  -y„R 

Y-*  cos  *e  1  e  0 

_  _ I 


2»(a1  +  i»e1)R3 
♦  Y0R  -  YoR2nB) 


1  -  YjR  sin  * 


.  “Y,b  -Y0R  I 

Y°5  Id  *  Y,R  sin  «C*f 

2ir(o1  +  io>ej)RH  1  | 


+  YqR2)  -  sin2  *(15  +  15y0R  ♦  6y2R2  +1 


II 


-Y,h  -y0R 

?  Y  f°  .  12  *  r,R  sin  ♦  *  t„R 

2.(0,  ♦  iac,)!1  1  >  ° 


(J 


+  YjR  sin  *)  -  sin2  *(3  +  3y0R  +  YqR2 


p. 


r.w'gwraatwg'fti 


Dipole 

Type 


£1 


Table  8.  Subsurf ace- to -Air  Propagation  Formulas 


x  (3  +  3y0z  +  y5z2) 


YjP  COS  + 

,  -Yih  ~Y0Z 
sxn  <ie  1  e  u 

2^(0^ 

+  iwejjz2 

x  (1  + 

Y0Z) 

Yfm  cos  + 

J  -Yih  ”Y0Z 
sxn  ipe  1  e  u 

2Jf(o^ 

+  iue^z2 

x  (1  + 

Y02> 

Y^m  sin  ifi 

cos  *e-Y>he-,»Z 

2ir(o 

+  ^l)*3 

x  (3  ♦  3y0z 

♦  YqZ2) 

•  .  •  .  -Yxh  -Y0z 

Yxp  sin  $  sin  ^e  1  e  0 

2*(Oj 

<•  iuej)z2 

X  (1  ♦  Y0z) 

- r 

9  •  •  “Ylh  “YnZ 

Y^m  sin  $  sin  ipe  1  e  u 

-Y,h  -yqz 

YjP  cos  $  cos  ^e  1  e  u 


2it(,al  *  iuejjz2 


*  (1  ♦  Y„z) 


,  -Y,h 

Y^m  cos  <{>  cos  ipe  1 


2ir(Oj  +  iue1)z 
x  (1  +  y„z) 


ft 


% 


■  *^fcb£  sagrawggg 


Formulas  When  z  »  p  (|n2j  >_  10,  z  >_  106,  z  >_  Sh) 
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VT»Z 


cos  te'Y*he-T°: 


i 

p  sin  *e  ^e  Y°Z 

2vy^z 

*  (1  +  Y0z  + 

3 

YqZ2) 

m  sm  4>e  ■ 

Y*v 

2*z3 

x  (1  +  YqZ  + 

V  ' 

-Y,h  -ynz 
s  le  1  e  0 


2tty,z3 


x  (1  ♦  T0Z  +  Yjz2) 


VY°Z 


YV 


ITY^ 

X  (3  +  3y„z  +  Y^Z2) 


2irYiZ 


x  (3  +  3y„z  +  YqZ2) 


he-Y0z 


-Yih  -Y0z 

m  sin  $  sin  *  cos  <lie  1  e 


2irz3 


x  (1  +  y0z  +  YqZ2)  x  (3  +  3yqz  +  YqZ2) 
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J  "Yjz  -Y0D 
VjP  cos  tye  1  e  0 


2ir(Oj  +  iuejjD2  (1  +  Y°DA) 


~Yiz  ~Yod 
m  cos  me  1  e  u 


-77-  y;  .  -,.-[(!  ♦  YjD  sin  *)  { 
Znfoj  +  iuejJD^  1 

sin2  «(»(15  +  15y0D  +  6y2D2  +  y3D3)] 


A  -Y,'z  -Y0D 
3  cos  $e  1  e  u 

2Tr(0j  -ue^D3 


(1  -  y:D  sin  <ji  +  y0D  -  y|D2nB) 


A  'Yiz  "YnD 
y.a  cos  $e  1  e  u 

— - =—• (1  +  y.D  +  y2D2A) 

2*(a1  +  iuCjJD3  0  0 


-Y,z  -YnD 


sin  Ae  1  e  0 


2Tr(Oj  +  iuejJD3 


[2  +  YjD  sin  ij> 


♦  y0DCl  +  A  +  yxD  sin  <0  -  sin2  ^(3 


Yjin  sin  $e  Y*Ze  Y°° 


2TTCOJ  +  iuejJD3  +  Yo°(1  + 
sin2  <i(3  +  3y_D  +  y2D2)] 


f  *****  Hi’**'-*  *  ^  1  •'*  AjZiM&f/aXilt  fetasi-E# 


■  w  «iv*v  £  i  •-*>  — 5t*J \  j»;  *■  *a. 
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Magnetic-Field  Air-to-Subsurface  Propagation  Formul'J 


Dipole 

Type 

HP 

- 3 

H*  ! 

_ 3 

VHD 

0 

1 

1 

| 

-y  z  -yftD  l 

p  COS  ).  '  e  °_(1  ,  DA)  I 

2irD2  0  S 

1 

■4 

% 

VMD 

-YiZ  -y0D 

-  m  COS  j^-pTr- - CC1  +  YjD  sin  *)(3  +  3y0D  ♦  y2D2) 

-  Sin2  <p(15  +  15yqD  +  6yqD2  +  YqD3)] 

1 

I 

1 

0  I 

| 

A 

HED 

•  x  ”Yi z  -Y0D 

p  sm  4>e  1  e  0  , 

*- - - - = — 1 - [2  +  YiD  sin  ip 

2irYlD3  1 

+  y0D(1  +  A  +  YjD  sin  ip)  -  sin2  ip(3  +  3y0D  +  YqD2)] 

1 

,  -YiZ  -Y0D  | 

-  PCOS£1D3e  -(1  '  ^  ^  *  *  Y"1 

1 

re 

‘1 

i 

1 

END 

°  **»  2  *  -r„DCi  *  A) 

-  sin2  $(3  ♦  3yqD  +  YqD2)] 

1 

1 

X  _Ylz  ~yod  1 

-  °  C0S  %D»  6  -  (1  *  V  *  | 

1 

i 

§ 

2  "Y0D 

3 - — (1  -  YjD  sin  4>  +  y0D  -  Y^nB) 


lW 


1  +  Y„D  +  Yo°2a) 


me 

2tt(y? 


‘V  -v 


r:^{a  +  *iD  sin  <10  I  (9  +  9y0D  ♦  4Y2D2  +  Y^3) 


-  sin2  *(15  +  ISyqD  +  6YqD2  +  YqD3)]  -  sin2  *(75  +  75y0D 
+  33YqD2  +  6y3D3)  +  sin4  *(105  +  105yqD  +  45yJjD2  +  6Y3D3)J 


sin  *  cos  *e  YjZe  Y°D, 


2ir(Yf  -  Y§)04 


(1  +  YjD  sin  *) (3  ♦  3y„D  +  YqD2) 


sin2  *(15  +  15y„D  +  6y§D2  +  y$D3)] 


m  sin  *  cos  *e  YjZe 


2ifYiD4 


[(3  ♦  3YoD  ♦  YgD2) 


-  sin2  *(1S  ♦  15y0D  +  6y2D2  ♦  y^D3)] 
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qjngaw^n^W 
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.'«»AJ-^1,!SSI!MII*JI^i-,iw. 
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Table  11.  Air-to-Subsurface  Propagation  Formulas  Whenf 


Dipole 

Type 

EP 

E* 

| 

Ez  % 

.•  ■% 

> 

VED 

YlPe-YiVY°p 

-  - - : — r-i*1  +  t0pA) 

2ir(Oj  +  lue^p2 

0 

1 

Ji 

1 

pe  YlZe  Y°P  1 

2ir(o,  +  iue,)p3  1 

.3 

1 

x  (1  +  y„P  +  YqP2A)  I 

1 

l 

VMD 

0 

me‘YiZe"Y0P 

-  TT - * — TT^1  +  Yih> 

2ir  (o  j  +  lue^p*  1 

x  (3  +  3y0P  +  YqP2) 

i 

1 

o  i 

i 

*& 

1 

HED 

x  *Ylz  "YnP 
p  cos  $e  1  e  u 

x  -Ylz  "Y0p 

p  sxn  *e  1  e  0  fo  . 

| 

p  cos  $e  Yl~e  Y°P  1 

2*^  +  iwe^p3 

x  (1  -  Yxh  +  y0p  -  YqP2JIB) 

2ff(a1  +  iuEjJp3 

+  YqPC1  +  A  +  Yjh)] 

2ir(o1  +  iwe^p3  I 

I 

*  [(3  +  3Y0p)sin  \l>  -  Af0p  +  YqP2B]  | 

HMD 

A  -Yiz  *Y0P 

Yjm  cos  $e  1  e  0 

Yjm  sin  $e  Y*Ze  Y°P 

X  "V  "V 
m  cos  $e  1  e  u 

2*^  ♦  iwej)p3 

x  (1  +  y0P  ♦  YqP2A) 

2*^  +  iuE^Jp3 

x  [2  +  y0pU  +  A)] 

2ir(0j  +  iucjjp2  ] 

x  (Y§)(1  +  Y0PA)  J 

in  Formulas  When  p  »  h  (|n2|  :>  10,  p  >_  106,  p  >_  5z) 


4  (1  ♦  Yxh) 


x  (3  +  3y„p  +  YfjP2) 


sin  Ae  Y*Ze  Y°P 


2ny jP  * 


[2  +  Yjh 


Of0P  +  YqP2b]  +  Y0P(1  +  A  +  Yxh)] 


«‘YiVYop 


2^r— a  -  Y0PA) 


2jty,p 


x  (1  -  Yah  +  y0P  -  YQP2nB) 


-YiZ  -Y0P 

m  cos  £e  1  e  0 


*  [2  +  Y0p(1  +  A)]  x  (1  +  YgP  ♦  Y§P2A) 
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me^V^ 


-  f ;  2  — zTT^1  +  ^ih) 
2iiCy2  -  Yq)p5 

X  (9  +  9y0P  +  4YqP2  +  YqP3) 


sin  Ae 


-Y,z  -Y„P 


*  'r>h) 


X  (3  +  3y„P  +  YaP2) 


•  *  *YiZ  -Y0P 

m  sm  Ae  1  e  u 


2»YiP 


*  (3  +  3y0P  +  YqP2) 
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x  (3  ♦  3Yoh  +  y§h2) 


-Y,Z  -Ynh 

p  cos  $  sxn  <ie  1  e  u 
27th2 

x  (l  +  Y0h) 


mm 


2iry1h3 

x  (3  +  3Yoh  +  Yq1*2) 


x  _^iz  ~y 
a  cos  6e  1  e 


x  (1  ♦  Y0h  +  Y§h2) 


-Y,Z  -Ynh 

m  sin  6  cos  *e  1  e  u 


2itYjh4 

x  (12  +  12Yoh  +  5yqIi2  +  y^3) 


& 


Dipole 

Type 

EP 

E* 

Ez 

4 

1 

VED 

-Y0P 

Yjpe  0 

2^(0^  +  iwe^p2 

0 

pe-YOp 

2iritoe0p3 

i 

i 

-1 

1 

a 

ss 

x  (1  +  Y0PF) 

x  (1  +  y0P  +  YqP2F) 

i 

i 

iup0me  Y°P 

i 

i 

»e-Y°p  1 

VMD 

0 

2if  (y2  -  Yq)p4 

0 

2,tYip4  S 

x  (3  +  3y0P  +  YqP2) 

X  (3  ♦  3y| 

p  cos  ♦e  Y°P 

,  -YnP 

p  sxn  $e  u 

iup0p  cos  ♦e  u 

p  sin  ♦el 

HED 

2n(o1  +  iwGj^p3 

2ir(a1  +  iuejjp3 

2iry  ]P2 

2lfYlPl 

x  (1  +  Y0P  +  yJp2F) 

x  [2  +  YqP (i  +  F)3 

x  (1  •+  Y0PF) 

x  [2  +  y| 

£! 

1 

\ 

*3} 

iu)yQm  cos  $e  Y°P 

x  -YqP 
iuu0a  sxn  ‘ 

x  "TnP 

iuy0m  cos  $e  u 

% 

i 

m  sin  4et 

»4D 

2mriP3 

*  (1  ♦  y0P  +  YqP2F) 

2rr1p3 

x  [2  +  Y0?C1  ♦  F)] 

2irp2 

x  (1  +  y0pF) 

2xp  H 

I 

x  [2  +  r 

i 

§ 

'  f. 
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tagation  Formulas  for  p  >_  106  (|n2|  >_  10,  z  =  h  =  0+) 


*8* 

% 


•  x  "YqP 
sin  *e  0 


2*yxP3 

2irYxP3 

2ir(Yf  -  Y§)P4 

*  [2  +  Y0P(1  ♦  F)] 

x  (1  +  y0P  +  YqP2F) 

x  (3  +  3yqP  +  YqP2) 

a  sin  ♦e’Y°P 

x  ‘Y0P 

b  cos  *e  0 

-Y0P 

m  sm  fe  0 

2*p3 

2wp3 

2ipf1P‘t 

*  [2  +  Y0P(1  ♦  F)3 

V  > 

*  (1  ♦  Y3P  +  YqP2F) 

x  (3  ♦  3yqP  ♦  YgP2) 
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